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Several arguments suggest parity violation may be observable in high energy strong interactions. 
We introduce new analysis tools to describe the azimuthal dependence of multi-particle distribu- 
tions, or "azimuthal flow." Analysis uses the representations of the orthogonal group 0(2) and 
dihedral groups Dn necessary to define parity completely in two dimensions. Classification finds 
that collective angles used in event-by-event statistics represent inequivalent tensor observables that 
cannot generally be represented by a single "reaction plane." Many new parity-violating observables 
exist that have never been measured, while many parity-conserving observables formerly lumped 
together are now distinguished. We use the concept of "event shape sorting" to suggest separating 
right- and left-handed events, and we discuss the effects of transverse and longitudinal spin. The 
analysis tools are statistically robust, and can be applied equally to low or high multiplicity events 
at the Tevatron, RHIC or RHIC Spin, and the LHC. 
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I. MANY FORMS OF STRONG PARITY VIOLATION 

Why should strong interactions conserve parity? Large violations of parity symmetry in strong interactions are 
possible by several independent means. Here we discuss new model-independent methods to probe strong parity 
violation at hadron colliders, including the Tevatron, LHC and RHLC. 



A. Strong Parity Violation at Low Energy 

Instanton approximations make one road to strong CP violations [1]. The model is described by a term 
QtiivufitriF 11 " ' F a P) added to the QCD Lagrangian. Here F Q/3 is the gluon field strength tensor, tr indicates the 
trace over the colors, and 9 is a parameter. Low-energy physics puts severe limits on the maximum size of such 
effects. Yet high energy collisions may well trigger a different vacuum phase with CP violation of order unity. Morely 
and Schmidt [2] long ago suggested testing strong P, CP, and T violations in heavy ion collisions. Recently the topic 
has seen great activity [3-10], including significant experiments we will discuss. Finch et. al.[ll] discuss some of the 
history of this development. 
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B. High Energy Strong Parity Violation 

High energies offer new possibilities for strong parity violation. As such, they are natural targets of RHIC and 
LHC experiments. Consider for example a higher-derivative model given by the QCD Lagrangian plus L e , where 

U = ^tr(e^F^F^F^). (1) 

In vector notation the Lq Lagrangian goes like c a b c E a x (B^ ■ B c ), where E a and B a are the color electric and magnetic 
fields, and c a b c are color group structure constants The Lagrangian represents a chromo-magnetic effect with three to 
six gluon vertices. Its transformation properties are odd under time-reversal, odd under parity and odd under CP. 
Very little is known to limit this Lagrangian. 

In contrast, 't Hooft's low-energy model goes like E a ■ B a . It is a total 4-divergencc naively dropping out of the 
integrated action, and giving "zero" observable effect on equations of motion. Its induced effects, expected from 
non-perturbative global subtleties of gauge fields, cannot really be estimated in perturbation theory. In comparison, 
the effects of Eq. 1 are straightforward signals that naturally grow with energy. The effects amount to parity- violating 
multi-gluon vertices that have every reason to affect high energy multi-particle production. 

The coupling of L 6 involves a mass scale M, which new physics considerations would estimate in the TeV range. 
In the regime of small momentum transfers of O(MeV), the ratio (MeV/TeV) 2 <~ 1CP 12 would greatly suppress its 
effects. It follows that constraints on the parameter M from low energy phenomenology such as the neutrons electric 
dipolc moment [12] are relatively weak for most beyond the Standard Model theories which can induce the operator 
in Eq. 1. Meanwhile the proposal of transient, local P or CP- violating phase not probed by low-energy observables 
may be even more natural when applied in a higher derivative context. High energies are generally needed to probe 
higher derivative terms and any new phases they may produce. No further motivation for the L 6 model is needed 
other than an effective theory with gauge invariance and symmetry under the rest of the Lorentz group. 

Given that many other models of high energy strong parity violation exist, it is surprising that few if any have been 
tested. A common assumption that strong parity violating effects should be "small" has little experimental support 
at collider energies, except by extrapolation. Testing strong P and CP symmetry at the new energy frontier is now 
a priority. 

C. Focus on Signals 

While models have their purpose, it would be premature to concentrate on models when testing parity symmetry. 
The work here is concerned with developing signals of strong parity violation that do not depend on models. We 
focus on the symmetry properties of observables that might discover parity violation from any source. Recently 
the STAR[13] experiment observed large effects in multi-pion azimuthal flow that are not consistent with certain 
parity-conserving simulations. The effects are described as being "a signal consistent with several of the theoretical 
expectations" sensitive to charge separation predicted by strong CP violation. But what did STAR really measure? 
The experiment measured a parity-conserving statistic. 

We will review how a parity-even statistic became tied to concepts of "odd parity fluctuations," sometimes leading 
to confusion over basic principles. We believe that reference to indirect, mo del- dependent deviations from a Monte 
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Carlo simulation cannot possible compete with direct confrontation of data with observables that absolutely preserve 
or violate a symmetry. As we will review in Section IIIC, the analysis long ago became embedded in an overly 
restrictive model - the "standard event plane formalism." That step unknowingly made parity-violating observables 
much harder to construct. It also greatly limited exploration of the full richness in multi-particle data with or without 
a focus on parity. In the model- independent framework we present, the scope and power of event-by-event and 
azimuthal flow statistics is not only retained, but extended to many new observables, including a new concept we call 
"shape sorting" (Section III D). 

Three-dimensional tests of parity symmetry tend to be complicated, and sometimes hinge on idealized three- 
dimensional detector symmetry. Here we focus on two-dimensional tests because they are simple, and because the 
naturally high statistical power of multi-particle observables makes many independent tests possible. Which two 
dimensional projection one should use is in principle a matter of choice and experimental details. Due to recent 
interest in parity violation tests using azimuthal correlations, here we will focus on observables constructed in the 
transverse plane. 

It turns out that defining and detecting parity violation in two-dimensional azimuthal flow is quite subtle. First, 
we will show that any quantity that is parity-odd by two dimensional tests is odd under ordinary three-dimensional 
parity. A proper two dimensional parity test procedure is sufficient to rule out parity symmetry, but not necessary: 
there also remains a number of independent three-dimensional tests. 

Any naive assumption that a single test would prove parity conservation is not true, due to deep geometrical 
features of the two-dimensional space. "Parity" in two dimensions is a rarely discussed subject that is more rich 
and complicated than the three dimensional rendition. By analyzing the transformation properties under 0(2), the 
orthogonal group in two dimensions, we not only discover new observable features of parity symmetry, but we also 
revise the parity- conserving measures of azimuthal flow. 

D. Outline of the Paper 

Besides testing parity with two dimensional quantities, the paper is concerned with the consistent two-dimensional 
statistical description of azimuthal flow. We find that inherent transformation properties cause significant revisions 
to standard formalism describing multi-particle data, whether or not parity symmetry is under discussion. 

The first step is group classification. The easy mapping of azimuthal moments into complex numbers produces 
distinct tensor transformation properties that have been hidden in plain sight. Transformation properties are crucial, 
because the distinction of different tensors is observable. Transformation under the dihedral group appears here, 
because parity in two dimensions needs it. 

The next step develops event-by-event statistics in general terms of conditional probabilities, also done in Section 
III. It is very important not to impose models that bias the experimental description. Section IV shows how to build 
variables capable of testing parity symmetry. Rather than singling out any "best" variable we show how to develop 
many consistent variables in order to let experimentalists determine what is "best." This Section also explains the 
classification of the parity-motivated statistics used by STAR. Section V discusses spin-dependent tests. 
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II. MANY REASONS 2D PARITY IS SO INTERESTING 

Details of collider experiments and detector geometries often make it difficult to test for parity violation using 
three-dimensional observables. It is thus more convenient to appeal to observables which "live" in two dimensional 
subspaces. In doing so, it is important to note that testing parity violation using two-dimensional observables should 
be consistent with three-dimensional parity violation. 

In three dimensions parity is represented by inversion, P 3 u = — l 3 x3- It is not necessary to invert all three axes to 
test parity symmetry, as we now show. 

Consider any operation D e 0(3) with determinant det(D) = —1. The operation of swapping x and y axes, 
with z fixed, is an example that can loosely be called "a parity transformation on a subspace." If symmetry of the 
Hamiltonian H under D fails, then [H, D] ^ 0. 

By construction one can always write D = P 3 dR, where R € 50(3). When D-symmetry fails, then 

[H, D] = [H, P 3D R]=P 3D [H, R] + [H, P 3D ]R^0, 

and then 

[H,P 3D }^0, (2) 

given rotational invariance [H, R] = 0. It follows that P3£>-symmetry is violated by finding any single example of 
parity violation on a subspace. 

Next: the map from three-dimensional data to two dimensions is not always trivial. Given a 3-momentum vector 
p = (p x , p y , p z ), remove the z component by projection: 

(Px, Py, Pz) -> (Px, Py)- 

Apply the same projection to the 3-dimensional parity operator P 3 d = —1.3x3: 

/ -1 \ 
■P 3 Z5 = 0-10 

V o o-i/ 

The 2x2 operator has determinant +1: it is a rotation, not "parity." The ordinary belief that parity is equivalent 
to an inversion fails! 

Parity is defined as the discrete subgroup with determinant det = — 1 of the orthogonal group. We seek a matrix P 2 d 
which represents two-dimensional parity and is rotationally invariant. Let M inv be a rotationally invariant matrix, 

a 13 



-1 
-1 



M; 

Calculate the rotated matrix 

R T (8) ■ M mv ■ R(6) = 



7 5 



cosO sinO \ I a (3 \ I cosO —sinO 

—sinO cosO I \ 7 5 I \ sinO cos6 

Here R(9) e 50(2) and superscript T stands for the transpose. Setting R T (0) ■ M inv ■ R(9) = M inv gives 6 = a, 
7 = —fi. All rotationally invariant matrices M inv then satisfy 

a (3 

M inv = | j (from rotational symmetry). (3) 

—j3 a 
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Compute the determinant: 



det(M inv ) = a 2 + 1 > 0. 



By requiring rotational symmetry, the condition det(P 2 D) = a 2 + (3 2 = — 1 is impossible. We find the remarkable fact 
that no real rotationally-invariant operator exists for 2D-parity. 

Candidates for 2D-parity must break rotational symmetry because all 2D matrices with det = — 1 correspond 
geometrically to a mirror-reflection about some axis in the plane. For example Po — diag(l, —1) reflects y-components 
around the x-axis. Any rotation of the matrix will also have det(R ■ Pq ■ R T ) — det(P ) = — 1 and define a different 
"parity." It follows that a continuous infinity of 2D-parity candidates exists: yet not one of them has the features 
expected in three dimensions. 

Despite these facts, it would be wrong to conclude that 2D parity does not exist as a consistent concept. The 
problem of "2-parity" is the lack of a rotationally invariant representation. But in deriving Eq. 2 we used rotational 
symmetry of the Hamiltonian, and not the same symmetry of the parity representation. There are many cases in 
physics where an operator is not invariant, but transforms in a known way. For example angular momentum is not a 
rotationally invariant quantity but we have rotationally invariant tests that it is conserved. For physical purposes of 
testing parity symmetry, we need invariant test criteria, whether or not the operator itself is invariant. That is, the 
representation of parity is a side issue compared to what we wish to test. 

Continuing: choose a candidate P 2 d € 0(2) by some axis convention. Let e be an operator to test D-parity 
symmetry; it will be odd under Pib'- 



Since P 2D Pid = 1 we find the trace tr(e) = 0. Impose rotational invariance of the test: R T ■ e ■ R = e. Rotational 
symmetry was solved by Eq. 3 for any invariant matrix M inv . tr{e) — implies parameter a — 0. Thus an invariant 
test exists. It is unique up to a scale (3 we set to unity: 



This is the 2-dimensional Levi-Civita symbol. It is the generator of SO(2). Since any generator is invariant under its 
own transformations, its rotational invariance is obvious. 

Continuing, the result did not depend on the choice of P-2D- Therefore e is odd under every det = — 1 transformation 
in 0(2), regardless of the reflection axis. Classifying the numerous representations is the topic of the next Section, 
using the dihedral group. The Section will show how our studies beginning with "parity" discover new features of 
azimuthal flow with or without parity symmetry. 



Classification of transformation properties is fundamental to "azimuthal flow." It is conventionally described as 
taking the Fourier transform of particle distributions in the azimuthal angle <j) defined relative to the z axis. By 
attending to transformation properties, we find there is much more in properly describing azimuthal flow. Extra 
attention to the group representations of 0(2) and setting up event-by-event statistics pays off by revealing new 
parity violating observables. 




III. AZIMUTHAL FLOW RECONSIDERED 
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A. Dihedrals 

The group 0(2) has infinitely many distinct discrete subgroups, called the dihedral groups of order N, denoted Dm- 
Each D N group is developed from two generators 1Z, V with the defining properties 

TZ N = 1, V 2 = l, vnv=n- 1 . 

The geometrical interpretation of 1Z is a rotation by 2ir/N . The interpretation of V is a reflection about an unspecified 
fiducial axis in the plane. The property VVSP = TZr 1 means the sense of rotations is reversed after the reflection. 
Although the generator algebra refers to no particular axis, it cannot be realized by matrices without selecting a 
preferred axis. 

Each group Dn has 2N elements (-Rjvfc, PNk), f° r k — 0...N—1, with det{R.Nk) — 1, det{PNk) = — 1, corresponding 
to certain rotations and reflections. For N = 3 an appropriate set is 

°)" R3I= (-~i-1)' R32= (?"-t) ; 

:)■-(??)■ M-l?) 

Wigner[f4] cites these matrices in his textbook as his first pedagogical example of an "abstract group," and later 
develops it as a homework problem, with no mention of "parity." 

The dihedral group is non-Abelian, except for the case of N — 1 , and group multiplication does not close until all 
elements are included. Because of that, one cannot make a consistent classification of two-dimensional "parity" by 
arbitrarily choosing one element. Nevertheless, by the previous relation to parity on 3-dimensions, any single violation 
of dihedral parity (D —parity) symmetry suffices to prove that 3-dimensional parity symmetry has been violated. As 
consistent, the invariant test Pjvfe ■ £ ■ PNk = —c holds for all possible values of N, k. 




B. Collective Angles 

Let $ be a collective angle obtained from an event with N particles in the final state. It is a random variable whose 
definition remains to be specified: it is not the same concept as the "reaction plane." The joint distribution of any 
pion (say) at angle <fr and * among a set of 3-momenta p\...p n is defined by 

dN f ,, ,, ,, dN 



d 3 pid 3 p2...d 3 p n 

X 8{4> - <j){p U p 2 ...p„)) <5(* - *(pi, P2-Pn))- 

The expression leaves open the definition of ^ = ^ (pi , p 2 . . .p n )) because there are many possibilities. The point 
is important, and developed below. Supposing angles are constructed properly, any distribution with rotational 
symmetry will have a trivial marginal distribution f(^f). Rotational symmetry then predicts a simple conditional 
distribution f(cj> | 

/(0l*) = ^^^.W-*)- 
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That implies the double-distribution does not need to be measured. It suffices to compute * event by event, and 
collect the <f> — ^ data into a single distribution. In effect events are "rotated to add" coherently. 
Such a distribution is conveniently expanded in a Fourier series: 



The angular momentum eigenstates are not eigenstates of dihedral parity, and vice-versa. That is consistent and 
necessary from operators that do not commute. Note index m of D m is matched to the Fourier mode. Conversely, 
one can expand in a basis of good D-parity eigenstates: 



Then parameters v m are even and a m are odd under under all elements P m k ■ With the understanding that the absolute 
value of angular momentum J is implied, we indicate good D-parity quantum numbers by | J\ p = J p . 

There is an instructive but wrong way to over-simplify the classification. An even (odd) transformation property 
of cosine {sine) functions follows from the simplistic transformation <p —<P- Yet that operation, which selects one 
particular origin of (p = 0, is highly coordinate-dependent. Indeed the whole distinction between sincj) and cos<p 
hinges on a coordinate origin for = 0, which is arbitrary. Basing decisions about "parity" on sine-versus cosine 
origin conventions cannot be physically meaningful. The over-simplification implicitly works under the assumption 
that only one mirror-plane about a reference point = defines "parity." 

It is somewhat more meaningful to examine the rotationally invariant functions of ~ ^- We have seen that 
D-parity involves reflection through a reference axis in the plane. Given a single pion at angle 0, and an event with a 
single collective angle ^F, it is natural to reflect through an axis defined by one or the other. But a single-reflection test 
can only confirm one (1) particular D-parity symmetry. All the other D-parity transformations not even examined 
might contradict it. For example reflection symmetry of the distribution across an axis 90° to ^ might well be violated 
even if reflecting across ^ itself was a good symmetry. 

Thus more work is needed. Let 0* be a reference point on the unit circle. Reflection of and * about the 
corresponding axis defined by 0* is 



/(0l*) = £A» e '' 



m 





in 



4> — > (j) = 4> + 2(</>» —(/))= 2(f>* — (f>; 



* -> = * + 2(0, - *) = 20, - *; 



1 Note the transformation is anti-unitary, which is a convention inherited from the complcxification of /, 
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Now sin(4> — \&) is correctly odd under all such reflections, regardless of axis orientations. That means it must be 
related to the operator e we found was unique. The relation is very simple. Define multiplets 

X0 = (cosfa sinfa); 
X\j/ = (cos^ , sin^f): 
then x^, ■ e ■ Xcj, = sin((f> — v &). 

As we discussed earlier, under P^k this transforms like 

X\fy • £ • X(p — Y .Ti|f • P^ k ■ £ ■ PjV/c ' Xrfc = —Xty ■ £ ■ Xfy. 

Once more the choice of Pjvfc reflection does not matter. 



C. Many Distinct Collective Angles 



Evcnt-by-event analysis is a breakthrough [f 5] whose repercussions are still being developed. Our work focuses on 
transformation properties, and does not depend on the details of how the event sample is selected. We assume there 
are cuts on the magnitude of transverse momentum, rapidity, particle species, etc. without discussion. We turn to 
reconsidering collective flow angles, so far represented by a single symbol 

The arithmetic mean angle (fa) is coordinate dependent, and transforms poorly. Standard procedures of "circular 
statistics" [f 6] map data into SO{2) representations using 

const (cos m , sinm^> m ) =^ (cos m fa, sin m fa). (4) 

i 

Here fa are angles of a data sample. In event-by-event analysis, the index i runs over a set of particle labels selected 
from each event, producing a number f m from each event. 

By construction each ty m parameterizes an SO(2) tensor of rank-m. As a consequence, each multiplet transforms 
like to units of 5*0(2) angular momentum. Upon rotating the </>-coordinate system by angle S, this means 



cosm^ m , \ I cosm^' m \ I cosmS —sinmS \ I cosm,^ r , 
sinm^ m I \ sinm^' m I \ sinmS cosmS J \ sinm^ r . 

Since they transform differently each collective angle ^ m is independent, and each represents new and distinct 
information about the underlying system. It is not correct mathematics to propose <I>i = *2 = ^z — as a "symmetry" 
of a physical system. 

Similarly, in three dimensions different tensors of rank I, 2, 3... cannot be identified as equivalent. Such tensors will 
be found event by event by multiplying momentum components pjPkPi---, taking moments, and making conditional 
distributions. In no way would we imagine that the tensor of rank-3 (say) was kinematically determined by the tensor 
of rank-2 or rank-4. 

Closely related, and coming from Eq. 4, each different ^ m has the defining property 

TO 
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The symbol = here means "identical by definition." The angle ^2 labels an axis, which transforms like a "stick" 
without oriented ends; by definition ^2 = ^2 + 7i\ An angle identical to itself in the class m = 3 does not represent 
anything geometrically like a vector, but something like the orientation of an ideal Mercedes-Benz symbol. 

The nearest equivalent to such properties in three dimensions comes from eigenvectors e a of real symmetric matrices, 
which have no sign. A pitfall of ordinary notation allows one to write expressions such as e a ■ v and believe it means 
something invariant. Each number e a ■ v is a scalar, yet ill-defined because e a ■ v = —e a ■ v when e a = — e a . The sign 
issue seldom causes problems, because it cancels out in expansions of the form v = e a (e Q • v). 

Needless to say, one is not going to find the sign of an eigenvector by doing a physics experiment. The lesson is 
that conventional notation may fail to incorporate the transformation properties symbols have from their definitions. 
We cannot find a record of this being noticed for collective angles. There are many consequences. In our notation 
symbol ty 2 is the angle for elliptic flow. Consider the expressions u sin((j) — ^2)", u cos((f> — ^2)", which appear to be 
well defined. They appear to be rotationally invariant under the shift of origin (j) — > </> + S, 'J 2 — > ^2 + S. Yet neither 
exists as a physical observable, because 

sin((j) — ^2) = sin(4> — ^2 — = —sin(<j) — \I/ 2 ); 
cos{<p — ^2) = cos{4> — ^2 — 7r) = —cos((f> — ^2)- 

It follows that an experiment seeking these moments will either get zero, or obtain a bias created by an improper 
usage of 

It may seem a paradox that related quantities such as "sm 2 (0 — ^2)" can be a good observable. Use the identity 

sin 2 (<t> - * 2 ) = ^ '-. (5) 

This is manifestly a function of 2* 2 , and 2^2 = 2^2 + 2tt is a safe periodic variable. Its relation to a second-rank 
tensor nature comes using index notation: 

sin 2 '(cp - #2) = (£*■£■ x^) 2 = xix^xlx 1 ^ e u . 

The product of four x makes a tensor ranging up to rank-4 which has been projected into a true scalar. 

These facts clarify puzzles coming from the standard event plane formalism [17-20]. It is a semi- classical model in 
which a single unobservable event plane angle ^rp is used for all 

Irp{4>) = ^2 v m cos((j) - * flP ) + a m sin(4> - Vrp). (6) 

m 

Here symbol = mean a defining fact of the model. The model can be implemented with classical simulations in 
which reaction planes and multi-particle collisions are generated numerically. The authors of Ref. [15, 21] recommend 
interpreting the different * m as "estimators" of fjjp, with discrepancies to be explained by fluctuations. That step 
defines their model. For model-independent analysis we recommend just the opposite, and exploring each ^ m as 
the independent tensor variable it is. Whether or not concerned with parity symmetry, eachangular tensor variable 
represented by "JmOught to be of physical interest, because each represents independent physical information. 



D. Event-Shape Sorting 

If one believed all collective angles *f> m were equivalent, then all event angular shapes would be boring, and more 
or less summarized by "elliptic flow." Let us shift emphasis to the correlations found in multi-particle distributions. 
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Correlations in the distribution of different collective angles describe meaningful "event shapes." 

Here is a simple example. An event's "splash pattern" can be characterized by first choosing a reference angle such 
as the orientation of elliptic flow. By construction it is D-parity-symmetric relative to itself. Each subsequent 
m th moment (a Fourier shape) has a particular offset relative to the second (elliptic) shape given by & m — ^2- The 

(2) 

magnitude ai, is a measure of the parity-odd strength of each particular mode. The sign is meaningful, and represents 
a sense of twist of ^ m relative to \I>2- 

0.2 
0.1 

-0.1 

-0.2 

Re(/ 2 ) 

FIG. 1. Real and imaginary parts of second Fourier moments sorted into Type-Minus (red online) and Type- Plus(blue online): 
see the text. Third moments were fixed. 

Rotational symmetry constrains the physics much less than one might believe. A general 3-variable distribution 
takes the form 

f(cf>, %, *,) =J2 /m, mi!m2 *(m+pmi+«)e lm * +,pmi ^ +I9m2 *' 

m 

This is more restrictive and more interesting than the naive periodic / = f((f> — cf> — ^2)- It also depends on 
more variables than the model of Eq. 6 is capable of emulating in principle. Thus for both parity-violating and 
parity-conserving circumstances, there will be advantages to using a completely model-independent formalism. 

To explore such possibilities, we propose systematically recording event-by-event statistics to make distributions 
conditional on more than one collective angle, which we call "event-shape sorting." We anticipate that strong in- 
teractions will develop interesting event shapes that grossly violate ordinary expectations. In particular, almost all 
individual events must violate D-parity. 

This may be surprising, because a naive interpretation of symmetry in quantum mechanics would expect that 
every single event be symmetric. However that is not observed. Instead, parity and rotational symmetry in quantum 
mechanics are generally realized by events that actually break the macroscopic symmetry badly. Events generally 
come in right- and left-handed versions. Parity symmetry, if true, predicts equal numbers of left- and right-handed 
cases. After sorting, the physics of azimuthal flow makes particular shapes. Exposing different cases for study requires 
sorting the data into parity-determined types before making distributions. 
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FIG. 2. "Event-shape sorting." Statistics of opposite D-parity are used to sort azimuthal angle data into classes labeled 
Type-Plus and Type-Minus. Sorted events make two distinct distributions, shown in polar plot as solid and dashed, where the 
radius at each angle equals the distribution. Data comes from Fig. 1 as described in the text. 

1. Simulation Example 

A Monte Carlo simulation illustrates the procedure. A large number (5000) of 4-mode event distributions (modes 
m = 0...3) was generated with flat random Fourier coefficients in the range keeping the distribution positive. Distribu- 
tions were generated so that the third moments were fixed. The second moments were then sorted into two classes: 
Type-Plus with Jm(/ 2 ) < and Im(f 2 ) < Re(f 2 ), (3787 events) and Type-Minus (3826 events) with Im(f 2 ) > and 
Im{f2) > -Re(/2). Fig. 1 shows a scatter plot of the Type-Plus and Type-Minus moments. This particular selection 
merely illustrates one possibility. Fig. 2 shows visually the Type-Plus and Type-Minus distributions accumulated 
over all events using a polar plot, and normalized by the number of events. The shape information represents average 
"splash patterns" coming out of the correlations in the events. Visual inspection shows the opposite _D-parity. The 
(unsorted) sum of these two distributions is symmetric about the x-axis: the sum would pass statistical tests of 
D-parity symmetry based on sin{4> — ^i) — > —sin(4> — "fi). 

Given the enormous multi-particle statistics of both the LHC and RHIC, we believe that many different combinations 
of event-shape observables will be statistically robust, and extremely interesting probes of strong interaction dynamics. 

IV. TESTING PARITY SYMMETRY 

We proceed to construct parity symmetry-testing observables systematically. Our focus is on general, model inde- 
pendent observables. 

The ordinary Fourier moments a m , defined relative to case by case, transform as odd-quantities under Z? m -parity. 
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Item 


Product 


Outcome 


Cartesian-Form 


Trig-Form 


Parity 


1 


1^ ® 1r 

A ^ B 


Oab 

AB 

2+ 

AB 


xa ■ &B 

G ^ X a X e> 

l/2(x^:r B + £^:r B ) - 1/2xa ■ x B 


cos{<j>A — 4>b) 
sin((f)A — <I>b) 
/ cos(4>a + 4>b) sin(4>A + 4>b) \ 

Tab = 

y sin{4>A + </>b) -cos{4> a + 4>b) J 


+ 
+ 


2 




0*. 

Oa, 
2a, 


e' 3 x\s 3 

X A ■ 8 

l/2{x\s 3 +x 3 A s t ) - 1/2xa ■§ 


sin(4>A - 4>s) 

COs(<j)A — <t>s) 


+ 
- 


3 


«ab ® 0± B 


U AB AB 


(XA ■ X B ) 2 , (xa ■ £ ■ X B ) 2 


cos(2(4>a - 4>b)) 


+ 


4 


«ab ® 0j B 


^AB AB 


(XA ■ X B )(XA ■ £ ■ Xb) 


sin{2(4>A - 4>b)) 




5 


^AB ® ^CD 


n+ 

U ABCD 


™ (, 1 AB 1 CD ) 


COs(<f>A + <pB ~ <f>C ~ 4>d) 


1 
) 


6 






tr\ 1 AB /* flP ) 


COS(q>A T <Pb &**rp) 


1 
] 


7 


2ab ® 2 Cs 


U AB Cs 


11 \ 1 A B 1 C s) 


COs((f>A ~\~ (f>B <j>C 0s) 





TABLE I. A survey of parity-even and -odd observables in two dimensions. 7a b , etc. are the traceless symmetric parts of the 
matrices x\x 3 B . The Item number identifies terms highlighted in the text. Combinations not listed are indicated by dots. 



The most simple such observable is the diagonal sine moment. However it vanishes event by event by definition: 

sinm {4n - $> m ) = cos(m$f m ) ^ sin(m(j>i) - sin(m^ m ) ^ cos(mcf> l ) = 0. (7) 

i i i 

Each Fourier mode- to already has to- fold rotational symmetry, as an eigenvector of R m k- Each is effectively rotated 
into cosine form by the rule defining 

The physics lies in the correlations between different modes. The "best" probe will be an experimental question. For 

(2) (2) 

dcfinitencss consider a model distribution / = f(2(j> — 2^2)- Let a m , % be the moments of the angular distribution 
relative to \E*2- The distribution is estimated by averaging over all data tfiij of all events J, numbering N to t- A short 
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calculation gives 

— — V sin m (<j> u - # 2 ) = \a { £> sin{ m (2* m - 2* 2 )); 

ij 

-!- J2 cosm(<t>u - * 2 ) - cos( m (2* m - 2* 2 )). (8) 

tot ij 

(2) 

Parity violation with a ^2 angular reference probe is signaled by a,„ 7^ 

Table I shows the classification of several quantities, which we now discuss. Recall our notation J p indicates 
the SO (2) angular momentum absolute value | J|, and D-Parity P, simply called "parity" here. The first moment 
xa = (cos(f)A, sirups) transforms like V^. Classify the direct products xa <£> &b- A short calculation shows they 
decompose under SO(2) into 

1 ± ® 1 ± -> 2+ +0+ +0". 
Notice this is quite different from the SO(3) rule 1(8)1 = 2 + 1 + 0. In more detail: 

• xa-xb = cos{4>a — 4>b) is the only even-parity invariant (0 + ) made from xa®%b- When xa — > % and xb — > xq, ± 
it has long been used to characterize directive flow: Item 1, Table 1. 

• xa ■ £ ■ xb — sin{<pA — 4>b) is the only odd-parity, 0~ available from 1^ (g> l|j. When (f>A, 4>b are the azimuthal 
angles of pions with the + and — charges, it is the odd-parity observable of Kharzeev and Pisarski[22] liberated 
from reference to a z-axis. (See Table 1, Item 2 and Section IV A.) 

• A spin 4-vector s^ of a Fermion projects into a 1 T multiplct s in two dimensions with D-parity intrinsically 
opposite to x. Then (Item 2, Table 1) x ■ s ~ 0~, and x ■ e ■ s <~ + . See Section V for more about transverse 
spin. 

• An interesting invariant (Item 6, Table 1) is made by reducing two spin-2 tensors to a new + invariant. For 
definiteness, make the J — 2 element from pions A, B, using the traceless symmetric tensor 

tab, « = _ m . m Sij . 

Construct a similar spin-2 tensor 7* from the collective flow multiplet x^. The angle ^ 2 describes the orientation 
of the tensor's principal eigenvector, mod^ir). The unique + from combining the two tensors uses the trace: 

2\ B ®2%^f 0+ B9 = tr(J~AB ■ T*), 

= cos{<j)A + 4>B-2^). (9) 

One important rule of thumb is that functions of vn\<f>— ™ 2 ^'m p only exist as invariant observables when mi = rri2 
and m2/m p — integer > 1. This is guaranteed when making invariants by contracting indices of Cartesian tensors. 

A. STAR'S Measurements 

Some time ago Kharzeev and Pisarski[22] (KP) proposed studying parity violation observing two final state particle 
pairs, such as ir + and ir~ with momenta k\, k 2 . The KP P 3 £>-odd statistic is z ■ (k\ x fc 2 ). 
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On general grounds, any Lorentz invariant cross section that violates parity symmetry must be 2 linear in the 4- 
dimensional symbol c^ap. If no polarizations are measured there must be four (4) independent 4- momenta, which 
can be taken as two beam momenta p±, p%, and two other detected momenta fcf , ■ The two beam directions are 
equivalent to an energy (timelike) direction = (E, 0, 0, 0) and a z-axis direction p*{ — p%. Then 

e^ aP p A p B kfk^ E(Pi - P*) ■ (h x fc 2 ). 

We call this the e- invariant; it is the Lorentz-invariant generalization of the KP invariant. As mentioned in the 
Introduction, the role of a z axis ("z" -symmetry) is part of verifying parity symmetry on three dimensions that does 
not enter the question of falsifying parity symmetry with two dimensional data. 

Kinematic constraints make it impossible to find the e- invariant in a 2 — > 2 perturbative sub-processes. When only 
four particles are involved, conservation of momentum allows any of the four momenta in Eq. 10 to be expressed as a 
linear combination of the other three. As a consequence, the contraction with the totally antisymmetric eV a P yields 
zero, a result we call "2 — > 2 epsilon no-go theorem." The theorem stands as a barrier to simplistic 2 — > 2 parton-level 
calculations that do not include initial or final state radiation, which carry away momentum, or more complicated 
final states. The theorem highlights a theoretical pitfall whereby Born-level calculations can yield null results, while 
inclusive reactions with many particles might do just the opposite. 

While important, the KP invariant is not unique, because any function of the momenta (or different momenta) 
can multiply it. Similarly, the most simple odd-D-parity quantity 0^ B = xa ■ £ ■ Xb, which is equivalent to the 
KP invariant, is not unique. It is a general fact that any invariant can be represented as a product of primitive 
invariants made from the smallest representation. When it comes to observables, however, the expectation of every 
moment-product is independent of all the rest: 

KbOab) ^ Kb)(Qab); 
(0£ B 0+ D <W * Kb)Kd){Obf), etc. (10) 

Knowing every moment of all products is equivalent to knowing the full distribution. It goes without saying that 
a model- independent approach makes no a priori relations between moments ("cumulant analysis"). Instead, each 
different moment contains information ready to explore. 

In Ref. [23] , Voloshin approached the question of parity-odd distributions of charged pions using the reaction-plane 
formalism. STAR[13] repeated the argument. The first step postulates pions correlated with the "y" transverse 
direction of angular momentum, described as 

^ = 1 + 2v 1 cos(Acf>) + 2v 2 cos(2A<f>) + 2a 1 sin(A(f>) + 2v 2 sin(2A(f>) + ... (11) 

Here A<f> = <f> — ^rp introduces the assumption of one universal reaction plane angle ^ rp. Refs. [13, 23] then turn 
to correlations of two particle species a, b of the form 

{cosA(j) a cosA(j) b ) - (sinA<p a sinA<j> b ) = (cos((f> a + (f> b - 2^ RP )) (12) 

-> (viv 2 ) - (a 2 ) + B in - B out . 



2 The statement assumes uppolarized beams 
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The last line is used to define B in — B out as any terms not coming from products of Eq. 11. The product of two 
sine functions led to interpreting (cos((f> a + (f>b — 2 1 $ i rp)) as a measure of charge-dependent parity violation. This 
introduces the concept of parity fluctuations. 

Experimentally STAR observed surprisingly large values of {cos((j) a + 4>b — 2^rp)). The STAR paper responsibly 
emphasizes that the correlation is actually even under parity, requiring no parity violation for its occurrence. Near 
the end of Ref . [23] , its authors also carefully noted parity-conserving resonance effects that could produce the same 
correlation. More recently, several papers [24-27] pointed out that azimuthal correlations stemming from transverse 
momentum and charge conservation can also explain the STAR data. In addition, effects of cluster particle correlations 
could also accommodate the data [28] . 

Nevertheless Ref. [23] maintains that the quantity is a measure of strong parity-violating fluctuations, adding 
emphasis that certain parity-conserving simulations do not reproduce it. STAR claims that "local parity violation 
cannot be significantly observed in a single event because of the statistical fluctuations in the large number of particles, 
which are not affected by the P- violating fields..." and "Improved theoretical calculations... are essential to understand 
whether or not the observed signal is due to local strong parity violation." 

We disagree. The perception that "local parity violation cannot be significantly observed" comes from the bias of 
the event-plane formalism and the focus on phenomenology of instanton approximations. Under that influence the 
role of parity fluctuations became central, despite their indirect and inconclusive nature. When exploring new physics, 
it is generally better not to rely on models, and return to the basic methodology of testing symmetries. In this sense, 
we find that testing strong parity violation is at least as important as testing strong parity fluctuation. 

The model-independent meaning of (cos(<p a + 4>b — ^2)) is given by the group classifications cited earlier, and in 
Table 1. There are no 0~ (parity odd quantities) involved. In coming from 2+ <g> 2^ (Table 1, item 6). We find the 
invariant interesting, and its measurement significant, but there is no model independent basis to claim it originates in 
parity violation or measures parity violation in any way. Note the observable is also not available from the invariants 
found in single-particle distributions. That may explain why simulation code tuned to single particle data did not 
reproduce its appearance. 

The question remains why the moment seemed related to parity violation in the first place. Trig-expanding of 
cos((j)a + 4>b — 2^flp) produces sin(<fi a — ^ Rp)sin{4>b — ^rp), among other terms. In the model of the RP formalism 
all tensors are labeled by a single ^rp, creating ambiguity. If ^rp — > "J 2 from elliptic flow, as used by STAR, then 
sin(A<fi a ) = sin(A(j) a + ir) = vanishes identically, and the same for b. To save the argument one might revert to 
^rp ~ ^lj which STAR mentions examining. Under angular momentum and parity, sin(<fi a — ^f\) 
where e is the 2x2 antisymmetric Levi-Civita matrix. This transforms like 0~, giving an appearance of parity-odd 
behavior, if standing alone. Yet by Eq. 7 it also vanishes at the single-particle, independence level where the argument 
is formulated. Recognizing this, some reference to fluctuations might be made, if only in a descriptive way we have 
not seen pinned down. Either way, the product of two such sines goes like X a ' £ ' X\$j Xfy ' £ ' X-fy . The product of two e 
obeys the identity 

£i 3 £k£ = SikSji - Sudjk- 

These delta functions exist from first principles, and are the "Clebsch" to project x a iXbjX^kX^/i — > 2+ b ® 2 J — > a b*. 
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V. 



SPIN DEPENDENT PARITY VIOLATING OBSERVABLES 



Transverse spin effects are large in high energy reactions, creating several puzzles. Studies of transverse spin effects 
are currently of great interest, and may well lead to new fundamental discoveries. 

Feynman's early treatment of transverse spin in deeply inelastic scattering [29] was flawed by an error in dealing with 
the quark mass. A well known and related mistake was introduced by Kane, Pumplin and Repko[30]. By examining 
low-order Feynman diagrams, the group concluded that any transverse spin effect in QCD must be proportional to the 
amplitude to flip a quark helicity, which is of order the quark mass divided by the energy m a /E in the diagram. The 
argument is wrong in general, and was never applicable in the conditions of small momentum transfers where applied, 
yet continues to be cited. Another argument based on Lorentz transformations was also popularized[31]. Under a 
boost along the z axis a Fermion spin 4-vector's transverse projection sV is invariant. The longitudinal projection 
s z transforms with a factor of the boost parameter 7, yielding the wrong conclusion that transverse spin effects have 
relative order m q /E once again. Meanwhile, it was known for decades [32] that transverse spin effects are inherently 
leading order in the parton model. The selection rules of chiral symmetry of hard perturbative interactions make it 
necessary to measure a chirally-odd distribution using another chirally-odd distribution [33]. The conditions of high 
energy reactions with several momentum transfer scales are much less demanding. There is actually no feature of 
QCD predicting any spin effects of high energy, soft strong interactions should be small. 

The STAR collaboration [34] measured large spin asymmetries of inclusive pp — > n° at y/s = 200 GeV. These 
measurements greatly extend the reach of the pioneering Fermilab experiment E704[35]. The size of these asymmetries, 
and even larger spin-correlations reported for 77 particles [36] shows once again that spin effect have no intrinsic high 
energy suppression. We concentrate here on collective flow analysis, noting that multi-particle reactions should not 
be thought equivalent to high-precision, hard-scattering collisions designed to test electroweak theories. 3 



Transverse spin greatly limits the number of possible new terms contributing to parity violating observables. Con- 
sider the single-spin case p\pb X, where p\ is the momentum of a polarized proton. By LSZ reduction, the proton 
initial state \pa, sa > is reduced to a Fermion interpolating field. Then the spin 4- vector sa enters the scattering 
amplitude only through the terms contracted with the on-shell spinor factors 



Thus only terms linear in sa are allowed. This severely constrains the possible parity-conserving and parity-violating 
terms that depend on spin. 

Turn to azimuthal flow, mapping by projection ^ Sa — (cos<p s , sincf) s ) in the transverse plane. Any general n- 
particle correlation with a single spin depends on some reducible tensor of the form 



A. Transverse Spin Classifications 



u(p A , s A )u(p A , s A ) = {f + m p ) 




2 



M i...jk 



(xix 2 ...x n ;s A ) = x\...x° n s k A . 



3 Electroweak-scale parity violation in longitudinally polarized proton collisions has been observed at PHOENIX [37, 38] and STAR [39] 
using the inclusive di-lcpton channel at ^/s ~ 500GeV 



17 



Consider a single-particle spin correlation for the purpose of illustration. Parameterize it with a tensor M. %3 , 

.VI" y.s',. 

M 13 decomposes into irreducible representations of 50(2) as 

1 / cos(tfi + 4> s ) sin{4> + 4> s ) 



sin(<j> + 4> s ) -cos{4> + </>g) 



+ ^cos(<p - <j> s )6 i3 ' - ^sin(<j> - <j> s )e i3 \ (13) 

where as before x — (cose/), sincj)). 

Since it must be linear in sa, an invariant distribution f(<j> \ </> a ) depends only on the invariant quantities associated 
with M 13 ' , which appear in the last line of Eq. 13. The distribution f((j)\(j) s ) of a single particle then must be of the 
form 

f((f> | 4> s ) = v + vicos ((p - (f) s ) + aisin {<p - <j> s ) . (14) 

Turn to parity. We first need to determine how three dimensional pseudo- vectors project onto two dimensions. 
For every 3D vector basis set {e a } one can define a pseudo-basis in the following way: 

Basis Pseudo — basis 



d ei x e 2 

e~2 e 2 x e~* 3 

e" 3 e" 3 x el 

Pseudo-vectors in three dimensions are linear combinations of the pseudo-basis. Define the vectorlike 2D basis by 
projection el — > el, el — >■ el, and the pseudovector-like 2D basis by projection: 

- ID . _ ID 

e\ — > e\ e 2 x e 3 — > e ■ e 2 

- 2D „ - 2D 

e 2 — > e 2 e3 x ei — 5- e • ei 

The pseudo- vector sa then projects into a pure pseudo-object sa as 

sa ^ sa = sie • el + s 2 e • e 2 (15) 
Use the fact that under dihedral parity the object e acquires a minus sign: P T e ■ P — > — e. Then 

cos(tf) — <p s ) — x ■ sa FNk > —cos{4> — <Ps) ~ 0~ 

sin(<f> — <p s ) — x ■ e ■ §a FNk > sin((f> — <f) s ) ~ + . (16) 

It is interesting this result is consistent with inversion and projection, i.e. three dimensional parity, while inversion 
itself is not "transitive" under projection. The analysis also shows a multiplicative rule for D-parity, with Eq. 13 
reading 1± <g> If -> 2~ s + 0~ s + 0+ . 

The presence of a non-zero coefficient associated with cos{<j)— <p s ) would indicate parity violation in a single-particle 
spin distribution. Numerous other examples can be developed. By the general rule that all invariants can be reduced 
to products of the most primitive ones, the D-parity violating moments take the form (0^ s 0jfc D ...0^ s ) . It is trivial 
to extend this to two spins. 
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B. Tensor Spin Correlation 

Finally, is the spin 2 irrep of Mij useful? It appears so. There are many ways of constructing scalar invariants from 
spin-2 objects. Consider the tensor T\ B for correlations of azimuthal directions of two final state particles (say 7r + 
and ir~ , Table 1, Item 1, bottom). Combine this with the spin-2 components of AAij, Eq. 13. There is one invariant 
we can construct out of the two spin 2 objects: 

2ab ® 2cs -> tr(T A B ■ Mcs) = cos(</>i + 4> 2 - <fe - <j> s ). (17) 

Despite the appearance of a cosine, this quantity is odd under D-parity. It is a parity-violating spin-dependent 
observable that could be measured by RHIC spin. 

Notice that it is impossible to construct an observable such as Item 6 in Table 1 from a single spin tensor correlation. 
The only possible combination is 

2 AB <g> 2 SS -> tr(TAB ■ M ss ) = cos(<t>i + <j) 2 - 2</> s ), 

but the tensor 2 SS is ruled out by the fact that only terms linear in s A are allowed. Other possible combinations 
can be determined using Item 5 in Table I, by replacing C with either A or B. A short calculation shows that the 
resulting observables are identical to items which can already be obtained from 1 ® 1 representations already listed 
in the table. 



C. Longitudinal Spin Classifications 

All our azimuthal observables can be extended to take into account longitudinal spin correlations. The general 
rule that spin must enter the distributions through a linear term only still holds, once again limiting the number of 
possible options. As before, the longitudinal spin parity violating observables can be constructed as 

Qab...n®xy...Zsi 

where A through Z can be any particles in the process. For instance, we can take any D-parity even observable 
involving N final state particles d + (<j>i, 4>2----4>n), and multiply it with the helicity of the initial state h(p A , s A ) to get 
a D-parity odd quantity: 

°12...JV L = d + (4> 1 ,4> 2 ,(j>N)h{p A , s A ) 

VI. SUMMARY AND CONCLUSIONS 

High energy parity violation is an exciting prospect. Until signals are observed in some form or other, and certainly 
when signals are observed and afterwards, it will be important to use an unbiased analysis framework that is free 
from model-based limitations. 

There are good reasons to seek high energy parity violation. When QCD is viewed within the modern framework 
of gauge-covariant derivative expansions, the parity conservation of the low energy gauge sector appears to be a 
kinematic artifact. That is because the ordinary QCD Lagrangian cannot violate parity without violating gauge 
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invariance. It is only by abandoning the most basic issues of boundary conditions in field theoriess that the 9- vacuum, 
9e.^vaptr{F^ F a @) term gets resuscitated for theoretical discussion. This is not always explained in the rush to 
construct motivation for the incredibly important task of testing parity symmetry. Since higher derivative effective 
actions do not have the same kinematic features, we find it very natural to look for parity violation at high energies, 
where effects are not necessarily small. Still, parity violation, like gold, will be where it is found, and nothing here 
favors one model over another. 

By systematizing the subtle facts of parity in two dimensions, we were led to the dihedral group, and from that 
to classification of statistics under the orthogonal group. A devious flaw of ordinary notation previously allowed 
distinctly different *f/ m to be identified with one another, concealing that each represents a distinct tensor, and each 
different f m = \t , m /(27rm). Attention to orthogonal group transformation properties makes this clear. In developing 
a model-independent description all ^ m are distinct a priori, and all moments are distinct a priori, with no conditions 
from "cumulants." 

Different \l/ m have been wrongly identified in the past as artifacts of statistical fluctuations. Thus the fascinating 
correlations of quantities defined with different f m remain for the most part unexplored. We have suggested a new 
concept of event shape sorting. It is the natural generalization of event-by-event statistics. Instead of imposing a 
model that all events are basically elliptical relative to one "unobservable reaction plane ^rp" , we suggest letting the 
data inform us of event shape categories. An interesting example concerns sorting events into right- and left-handed 
classes, which we illustrated by numerical simulation. 

The analysis tools are statistically robust because they involve nothing but moments of distributions. A huge 
amount of data is available from current colliders, suggesting that many independent studies can be conducted We 
anticipate that the classifications and analysis tools can be applied equally to low or high multiplicity events at the 
Tevatron, RHIC or RHICspin, and the LHC. 
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